Mapping the large-angle deviation from Gaussianity in simulated CMB maps 
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The detection of the type and level of primordial non-Gaussianity in the CMB data is essential to 
probe the physics of the early universe. Since one does not expect that a single statistical estimator 
can be sensitive to all possible forms of non-Gaussianity which may be present in the data, it is im- 
portant to employ different statistical indicators to study non-Gaussianity of CMB. In recent works 
we have proposed two new large-angle non-Gaussianity indicators based on skewness and kurtosis 
of patches of CMB sky sphere, and used them to find significant deviation from Gaussianity in fre- 
quency bands and foreground-reduced CMB maps. Simulated CMB maps with assigned type and 
amplitude of primordial non-Gaussianity are important tools to determine the strength, sensitivity 
and limitations of non-Gaussian estimators. Here we investigate whether and to what extent our 
non-Gaussian indicators have sensitivity to detect non-Gaussianity of local type, particularly with 
an amplitude within the 7 yr Wilkinson Microwave Anisotropy Probe (WMAP) bounds. We make a 
systematic study by employing our statistical tools to generate maps of skewness and kurtosis from 
several thousands of simulated maps equipped with non-Gaussianity of local type of various ampli- 
tudes. We show that our indicators can be used to detect large-angle local-type non-Gaussianity 
only for relatively large values of the nonlinear parameter f}$£ . Thus, our indicators do not have 
enough sensitivity to detect deviation from Gaussianity with the nonlinear parameter within the 7 
yr WMAP bounds. This result along with the outcomes of frequency bands and foreground-reduced 
analyses, suggests that non-Gaussianity captured in the previous works by our indicators is not of 
primordial origin, although it might have a primordial component. We have also made a compar- 
ative study of non-Gaussianity of simulated maps and of the full-sky WMAP foreground-reduced 
internal linear combination (ILC)-7 yr map. An outcome of this analysis is that the level of non- 
Gaussianity of the ILC-7 yr map is higher than that of the simulated maps for /{j ™ 1 within WMAP 
bounds. This provides quantitative indications on the suitability of the ILC-7 yr map as Gaussian 
reconstruction of the full-sky CMB. 

PACS numbers: 98.80.Es, 98.70.Vc, 98.80.-k 



I. INTRODUCTION 

The statistical properties of the temperature 
anisotropies of CMB radiation and of the large-scale 
structure of the Universe offer apowerful probe of the 
physics of the early universe [lJ-Q - Thus, for example, 
although a convincing detection of a significant level of 
primordial non-Gaussianity of local type (/nl^ 3> 1) in 
CMB data would not rule out all inflationary models, 
it would exclude the entire class of single scalar field 
models regardless of the form of kinetic term, potential, 
or initial vacuum state P, @. Such a detection, however, 
would be consistent with some alternative models of 
the physics of the early universe, and also with other 
nonstandard inflationary scenarios. On the other hand, a 
convincing null detection of primordial non-Gaussianity 
of local type (/nl^ — 0) w °uld rule out the current 
alternative models of the primordial universe (see, e.g., 
Refs. UflU). Thus, a detection or null detection of 
primordial non-Gaussianity of local type in the CMB 
data is crucial not only to discriminate or even exclude 
classes of inflationary models, but also to test alternative 
scenarios, offering therefore an important window into 
the physics of the early universe. 

However, there are various effects that can produce 
non-Gaussianity. Among them, the most significant 



are possibly unsubtracted diffuse foreground emission 
0, Hj, unresolved point sources @, possible systematic 
errors [l(|, and secondary anisotropies such as gravita- 
tional weak lensing and the Sunyaev-Zeldovich effect (see, 
e.g., Refs. [l|, [Ij])- 1 In this way, the accurate extraction 
of a possible primordial non-Gaussianity is a challenging 
observational and statistical enterprise. 

In the search for non-Gaussianity in CMB data, dif- 
ferent statistical tools can, in principle, provide informa- 
tion about distinct forms of non-Gaussianity. On the 
other hand, one does not expect that a single statisti- 
cal estimator can be sensitive to all possible forms of 
non-Gaussianity that may be present in CMB data. It 
is therefore important to test CMB data for deviations 
from Gaussianity by using different statistical tools in 
order to quantify or constrain the amount of any non- 



1 For most of these effects there seems to be no evidence of signif- 
icant non-Gaussian contamination within Wilkinson Microwave 
Anisotropy Probe (WMAP) sensitivity. Nevertheless, this is not 
guaranteed to hold true for the PLANCK experiment [ill ], due 
to its much higher sensitivity Q. Cosmic variance further com- 
plicates the problem, since some non-Gaussianity may arise from 
the uniqueness of the observed CMB sky. Deviation from Gaus- 
sianity may also have a cosmic topology origin Refs. [1^1 . 
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Gaussian signals in the data, and extract information on 
their possible origins. This has motivated a great deal 
of effort that has recently gone into the search for non- 
Gaussianity in CMB maps by employing several statis- 
tical estimators (an incomplete list of references is given 
the recent reviews [2-|j|, and also in Refs. [l3l - [l5j ) 

In recent papers [l3l IbH ] we have proposed two new 
large-angle non-Gaussianity estimators, which are based 
upon the skewness and kurtosis of the patches (spherical 
caps) of the CMB sky sphere. By scanning the CMB 
sphere with evenly distributed spherical caps, and calcu- 
lating the skewness S and kurtosis K for each cap, one 
can have measures of the departure from Gaussianity on 
large angular scales. Using this procedure we have car- 
ried out analyses of large-angle deviation from Gaussian- 
ity in both band and foreground-reduced maps with and 
without the KQ75 mask recommended by the Wilkinson 
Microwave Anisotropy Probe (WMAP) team for Gaus- 
sianity studies of CMB data [13l Ibil ]. Thus, we have 
found, for example, significant departure from Gaussian- 
ity in all full-sky band maps, while for the Q, V, W 
maps with the KQ75 mask, the CMB data arc consis- 
tent with Gaussianity. The K and Ka maps, however, 
show an important deviation from Gaussianity even with 
a KQ75 mask. For the full-sky foreground- reduced [l6r - 
[l9l | we have found a significant deviation from Gaussian- 
ity which varies with the cleaning processes, as measured 
by our indicators fl4j . 

Simulated CMB maps equipped with assigned primor- 
dial non-Gaussianity are essential tools to test the sensi- 
tivity and effectiveness of non-Gaussian indicators. They 
can also be used to study, e.g., the effects of foregrounds 
and other non-Gaussian contaminants, and also to dis- 
close potential systematics. In a recent paper a new al- 
gorithm for generating non-Gaussian CMB temperature 
maps with non-Gaussianity of the local type has been 
devised [2(|. In the simulated maps the level of non- 
Gaussianity is adjusted by the dimensionless parameter 
/^L al - A set of 1 000 CMB temperature simulated maps 
with the resolution of PLANCK mission for open 
(unfixed) values of f^™ 1 was ma de available [2(| . 

A pertinent question that arises is whether the indi- 
cators proposed in Ref. [llj, and used in Refs. [HI, 03] 
to detect non-Gaussianity in CMB data, have sufficient 
sensitivity to detect deviation from non-Gaussianity of 
local type with an amplitude f^ 31 within the bounds 
determined by the WMAP team in their latest data re- 
lease [2l|. Our primary aim in this paper is to address 
this question by extending the results, and by comple- 
menting the analyses of Refs. [HI, [lij in three different 
ways. First, instead of using CMB data, we use our sta- 
tistical indicators to carry out an analysis of Gaussianity 
of simulated maps equipped with non-Gaussianity of lo- 
cal type. Second, by using simulated maps with different 
amplitudes /},™ al , we make a quantitative analysis of the 
degree of sensitivity of our skewness and kurtosis indica- 
tors to detect primordial non-Gaussianity of local type. 
Third, we make a comparative study of the degrees of 



non-Gaussianity of simulated maps with different /nl s 
and the full-sky foreground-reduced 7 yr WMAP internal 
linear combination (ILC-7 yr) map [16j . An interesting 
outcome of this comparative analysis is that the level of 
non-Gaussianity of the ILC-7 yr is considerably higher 
than that of the simulated maps for f^l &1 within ob- 
servational bounds [2lj . This renders information about 
the suitability of the ILC-7 yr foreground-reduced map 
as Gaussian reconstruction of the CMB full-sky. 

II. NON-GAUSSIANITY AND SIMULATED 
MAPS 

A. Primordial non-Gaussianity of local type 

The first important ingredient in the study of non- 
Gaussianity is the primordial gravitational curvature per- 
turbations C(x, t), which were seeded by quantum fluc- 
tuations in the very early universe. In the linear or- 
der, the primordial curvature perturbation is related to 
Bardeen's curvature perturbation <&(x, t) [22j in the mat- 
ter dominated era by £ = (5/3)$ |23j |. The relation 
with CMB observations is given in the Sachs- Wolfe limit, 
where AT/T = - $/3 = -C/5 holds 0]. 

The lower order statistics able to distinguish non- 
Gaussian from Gaussian distributions is the three-point 
correlation function. Primordial non-Gaussianity can 
then be described in terms of the three-point correla- 
tion function of the curvature perturbations ^(x) or of 
its Fourier transform ^(k) by using the ensemble aver- 
age. Thus, the three-point correlation function counter- 
part in Fourier space — the so-called bispectrum — takes 
the form 

($(k 1 )$(k 2 )$(k 3 )) =S 3 (k 123 )B <i (k 1 ,k 2 ,k 3 ) , (1) 

where <5( 3 )(ki 23 ) = <5^ 3 ^(ki + k 2 +k 3 ) enforces that the 
wavevectors in Fourier space have to close to form a tri- 
angle, i.e., ki + k 2 + k 3 = 0. Thus, the bispectrum is a 
function of the triplet defined by the magnitude of the 
wave numbers (ki,k 2 , k 3 ). 

In the examination of primordial non-Gaussianity, the 
bispectrum of curvature is rewritten in the form 

B $ (k 1; k 2 ,k 3 ) = ($(ki)$(k a )$(kB)) 

= / NL (2tt) 3 <5 3 (k 123 ) F(k u k 2 ,k 3 ), (2) 

where / NL is an overall amplitude (dimensionless) pa- 
rameter, 2 which can be constrained by the CMB data, 
and where F(ki, fc 2 , £3) is the so-called shape of the bis- 
pectrum and encodes the functional dependence of the 



2 The subscript 'NL' stands for nonlinear. This notation arises 
because an often used phenomenological parametrization for the 
non-Gaussianity of 3>(x) can be written as a nonlinear transfor- 
mation of a Gaussian field (see below for more details). 
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bispectrum on the triangular configurations. The bis- 
pectrum shape is used to specify the configuration of 
the wave vectors of curvature perturbations that give the 
highest contributions to the amplitude of the bispectrum. 
In this way, the bispcctra are usually classified according 
to shapes of the triangles that give rise to the highest 
non-Gaussian signal. The most studied shapes are (i) 
local, where the bispectrum is maximum on the squeeze 
triangle configuration, k\ w &2 3> ^3! (ii) equilateral, 
with the bispectrum maximized on the equilateral limit 
ki m fc 2 ~ ^3! and (iii) orthogonal, where the bispectrum 
has a positive peak at the equilateral configuration and a 
negative valley along the flattened (elongated) triangular 
configuration k% « fci + ki [![. Predictions for both the 
amplitude / NL and the shape of B^(k\,k2,k^) depend 
on the early universe models, which makes apparent the 
power of the bispectrum as a tool to probe models of the 
primordial universe. 

In this paper we shall deal with non-Gaussianity of lo- 
cal type, which is the most studied type of deviation from 
Gaussianity. For this type of primordial non-Gaussianity 
the curvature perturbation in the real space can be split 
into two components: the linear term $l (representing 
the Gaussian component) plus a non-Gaussian second 
term, namely 3 

$(x) = $ L ( x) + ^° L cal ( *£(x) - ( *£(x) ) ) , (3) 

where both sides are evaluated at the same spatial loca- 
tion x. 

Finally, we recall that the latest WMAP constraints 
reported on /jSjL al at the 95% confidence level are given 
by -10 < /^°£ al < 74 [2l|. These bounds will be used in 
our analyses in one of the following sections. 



B. Simulated CMB maps 

The first simulated CMB temperature maps endowed 
with primordial non-Gaussianity introduced through a 
non-Gaussian parameter /nl were generated by the 
WMAP team and reported in Ref. [g. The WMAP 
algorithm was generalized so as to improve computa- 
tional speed and accuracy, and also to include polariza- 
tion maps in Ref. |26| fs ee also Ref. j27|). More recently, 
Eisner and Wandelt [2(| presented a new algorithm and 
generated 1 000 high-angular resolution simulated non- 
Gaussian CMB temperature and polarization maps with 
non-Gaussianities of the local type, for which the level of 
non-Gaussianity is determined by the parameter /^l ■ 
In their algorithm a simulated map with a desired level 
of non-Gaussianity f^ff 1 * s sucn that the spherical har- 



monic coefficients are given by 

_ L I /'local NL / a 

aim — a e m + /nl ' a im-> V±) 

where 

^ im and a £ nl are, respectively, the linear and 
nonlinear spherical harmonic coefficients of the simulated 
CMB temperature maps. 

In the next sections we shall use the set of 6, 000 CMB 
temperature simulated linear and nonlinear component 
maps, generated for an arbitrary (unfixed) value of /jjl, 
which was made available in Ref. [2(| 4 with the resolution 

of the PLANCK mission 

III. NON-GAUSSIANITY INDICATORS AND 
ASSOCIATED MAPS 

To make our paper clear and self-contained, in this 
section we describe the two statistical non-Gaussianity 
indicators and their associated maps, which can be cal- 
culated from either simulated or real CMB temperature 
(input) maps. The procedure delineated here will be used 
in the following sections to investigate large-angle devia- 
tion from Gaussianity. 

The most important underlying idea in the construc- 
tion of our non-Gaussianity indicators and the associated 
maps, is that one can access the deviation from Gaussian- 
ity of the CMB temperature fluctuations by calculating 
the skewness S and the kurtosis K, which measure, re- 
spectively, the symmetry about the mean temperature, 
and the non-Gaussian degree of peakness of the temper- 
ature distribution. 

Clearly the calculation of S and K from the whole 
CMB sky sphere of temperature fluctuations data would 
be a crude approach to measure the complexity of the 
non-Gaussianity of the CMB map. However, instead of 
having just two dimensionless numbers, one can go fur- 
ther by taking a discrete set of points j = 1, . . . ,N C ho- 
mogeneously distributed on the CMB sky sphere S 2 as 
the center of spherical caps of aperture 7 (say) , and cal- 
culate Sj and Kj for each cap j of the CMB temperature 
sky sphere. The values Sj and Kj can then be taken as 
measures of the non-Gaussianity in the direction (6j, 4>j) 
of the center of the spherical cap j . Such calculations for 
the individual caps thus provide quantitative information 
(2N C values) about non-Gaussianity of the CMB data. 

The above procedure is a constructive way of defin- 
ing on S 2 two discrete functions S(6, 4>) and K(9, <f) that 
measure departure from Gaussianity of a CMB temper- 
ature (input) map. In other words, these functions can 
be constructed from an input CMB map through the fol- 
lowing steps [l3| : 

i. Take a discrete finite set of points j = l,...,N c 
homogeneously distributed on the sky sphere of a 



3 It should be note, however, that this is not the only way to pro- 
duce the bispectrum of local type. Multi-field curvaton inflation, 
for example, can also produce bispectrum of local type l2- r l . 
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FIG. 1: The left panel shows a simulated Gaussian CMB map (/JjL al — 0), while the right panel depicts a non-Gaussian 
simulated CMB map calculated for /{j ™ 1 = 5 000. Temperatures are in mK. 



CMB input map as the centers of spherical caps of 
a chosen aperture 7; and calculate for each cap j 
the skewness and kurtosis given, respectively, by 

Np _ 

S ^F7E( T '- r ) 3 ' ( 5 ) 
iV P CT j i=i 

Np _ 

K i = ^iE( T *- T ) 4 - 3 ' (6) 

iV P CT i i=i 

where T, is the temperature at the i th pixel, 7} 
is the CMB mean temperature of the j th spheri- 
cal cap, Np is the number of pixels in the cap j , 

and a 2 = (1/N p ) J2 t Ii {T t ~T) is the standard 
deviation. 

Clearly, the whole set of values Sj and Kj (for all j) 
obtained through this discrete scanning calculation 
process, along with the angular coordinate of the 
center of the caps (9j , <f>j ) can be taken as measures 
of non-Gaussianity in the directions of the center of 
spherical caps j . 

ii. Patching together the Sj and Kj values for all 
spherical caps j, one has two discrete functions 
S = S(9, </>) and K = K(9, <t>) defined over a two- 
sphere S 2 . These functions provide measurements 
of the non-Gaussianity as a function of (9, 4>). The 
Mollweide projections of the functions S = S(9, (f>) 
and K = K(9, </>) are nothing but skewness and 
kurtosis maps (hereafter denoted by S map and K 
map). 

Clearly, the functions S = S{9, 0) and K = K{9, 0) 
can be expanded into their spherical harmonics. Thus, 
for example, for S = S(9, cf>) one has 

00 1 

S(0,^) = J2 E bt m Y tm (6,tt>) , (7) 

i=0 m=-l 

and the corresponding angular power spectrum 



which can be used to quantify the amplitude (level) and 
angular scale of the deviation from Gaussianity. In this 
paper we shall use the power spectra S( and Kg to esti- 
mate the departure from Gaussianity and calculate the 
statistical significance of such a deviation by compar- 
ison with the corresponding power spectra calculated 
from S and K maps obtained from input Gaussian maps 

(/NL Cal = 0). 

In the next section we will use these statistical indi- 
cators to carry out analyses of Gaussianity of simulated 
maps endowed with non-Gaussianity of local type of dif- 
ferent levels to test the sensitivity of S and K indicators 
to detect primordial non-Gaussianity of local type. Fur- 
thermore, we will make a comparative study of the de- 
grees of non-Gaussianity of different simulated maps and 
the ILC-7 yr WMAP map. 



IV. NON-GAUSSIANITY AND SENSITIVITY 
OF S AND K INDICATORS 

A. Analyses and results for simulated maps 

It is clear from the previous section that in order to 
study the sensitivity of the non-Gaussianity indicators S 
and K, that is, to construct the functions S — S(9, 4>) 
and K = K(9, </>) and the associated S and K maps, we 
need CMB input maps. The input maps used in our anal- 
yses in this section are high-angular resolution-simulated 
CMB temperature maps endowed with non-Gaussianities 
of the local type introduced through different values of 
the dimensionless amplitude parameter Jnl^'- 

Figure [TJ shows two examples of such simulated CMB 
maps. The left panel gives a Gaussian map /^L al = , 
while the right panel shows a non-Gaussian map for 
/k.£ al = 5 000 . We have taken these values for f^ al 
as an illustrative example that makes the non-Gaussian 
effects visible to the naked eye through the comparison 
of these simulated maps. 

Figure [5] gives an illustration of typical skewness S (left 
panel) and kurtosis K (right panel) maps, generated from 




FIG. 2: Examples of skewness (left panels) and kurtosis (right 
respectively, /}, L cal = 500 (first row) and /n™ 1 = (second row' 



input CMB simulated maps for /^L al = 500 (first row) 
and f^l al = (second row). The non-Gaussian maps 
show spots with higher and lower values of S(9, <fi) and 
K(Q,<j)), which suggest large-angle dominant multipole 
components (low €) in these maps. 5 In our calculations 
of all S maps and K maps, to minimize the statistical 
noise, we have scanned the CMB sphere with N c = 3 072 
spherical caps of aperture 7 = 90°, centered at points 
homogeneously distributed on the two-sphere. 

To obtain quantitative large-angle-scale information, 
as a first test of the sensitivity of S and K indica- 
tors, for each of the following values of the nonlinear 
parameter -/^°£ al = 0,500,1000,5 000- we have gener- 
ated 1000 simulated CMB maps, totalizing 4 000 sim- 
ulated maps. From each set of 1 000 simulated input 
CMB maps (fixed /nl^ f° r each set) we have calculated 
1 000 S maps and 1 000 K maps, from which we com- 
puted the associated power spectra, namely, S\ and K\, 
where i = 1, ••• ,1000 is an enumeration index, and 
I = 1, • • • , 10 is the range of multipoles we have fo- 
cused on in this paper. 6 The low £ multipole mean com- 
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panels) indicator maps calculated from simulated maps with, 



ponents Si and Kg arc then obtained by averaging over 
1 000 power spectra calculated from simulated maps with 
the different values of /nl 3,1 - The statistical significance 
of these power spectra is estimated by comparing the 
values of Si and Kg obtained from input maps generated 
for f^ff = 500, 1000,5 000 with the values of the cor- 
responding power spectra Si and Kg obtained from the 
Gaussian simulated map (/nl^ 1 = 0). 

Let us describe with some details the calculations of an 
average power spectrum. For the sake of brevity, we fo- 
cus on the skewness indicator S along with a set of 1 000 
simulated non-Gaussian input maps computed for a non- 
linear parameter fjS£ = 500, for example. 7 Following 
the two step procedure described in Sec. IIII1 from these 
1 000 simulated CMB maps we calculated 1 000 S maps, 
from which we computed 1 000 power spectra S^ in or- 
der to have the average values Si = (1/1000) J^I^a ■ 
From this MC process we have at the end ten mean mul- 
tipole values Si (£ = 1, •••,10), each of which is then 
used for a comparison with the corresponding average 
multipole value Si calculated by a similar procedure from 
Gaussian simulated maps (/nl^ = 0). This allows the 
evaluation of the statistical significance of Si by quanti- 



We have also calculated similar maps from the simulated maps 
for the other values of /jfjj al - These maps, however, provide only 
qualitative information, and to avoid repetition we only depict 
maps of Fig. [2] for illustrative purpose. 

The values of Imax for Si and Ki depend on the resolution of S 
and K maps, which clearly depend upon the number of spherical 



caps used in the scanning process. For N c = 3 072 one can go up 

to Imax = 45. 

7 A completely similar procedure can clearly be used for the kur- 
tosis indicator K along with other values of /^l^ ■ 
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fying the goodness of fit for Se with f^ al ^ and Se 
calculated for /^f^ 1 — (Gaussian maps). 

Figure [3] shows the average power spectra of the skew- 
ness Si (left panel) and kurtosis Ke (right panel), for 
i = 1, ■■• ,10, calculated from simulated Gaussian 
( /nl^ = 0) maps, and from CMB maps equipped with 
non-Gaussianity of the local type for which = 
500, 1000, 5 000. The 95% confidence level, obtained 
from the S and K maps calculated from the Gaussian 
CMB simulated maps, is indicated in this figure by the 
dotted line. 

To the extent that the average Si and Kg obtained from 
input simulated CMB maps endowed with /^f^ 1 = 500 
are within 95% Monte-Carlo (MC) average values of Se 
and K e for /^ al = 0, Fig.[3]shows that our indicators are 
not suitable to detect primordial non-Gaussianity of local 
type in CMB maps smaller than this level. However, this 
figure also shows that they can be effectively employed 
to detect higher levels of non-Gaussianity of local type. 
These results square with the statistical analysis we shall 
report in the remainder of this paper, particularly with 
the x 2 test of goodness of fit. 

To have an overall assessment power spectra Se and 
Ke, calculated from the input simulated non-Gaussian 
maps equipped with primordial non-Gaussianity with 
nonlinear parameter f^ 1 = 500, 1000, 5 000, we have 
made a x 2 test to find out the goodness of fit for Se 
and Ke multipole values as compared to the MC mean 
multipole values obtained from S and K maps of the 
Gaussian (/nl = 0) simulated maps. In each case, 
this gives a number that quantifies collectively the de- 
viation from Gaussianity. For the power spectra Se and 
Ke we found the values given in Table Q] for the ratio 
(reduced x 2 ) x 2 /dof (dof stands for degrees of freedom) 
for the power spectra calculated from non-Gaussianity 
of local type with = 500, 1000, 5 000. Moreover, 

the greater the reduced x 2 values (hereafter denoted sim- 
ply by x 2 )) the smaller the x 2 probabilities, that is, the 
probability that the multipole values Se and Ke and the 
mean multipole values obtained from the Gaussian maps 

(/nT 1 = °) a g ree - 



Non- Gaussian Parameter 


X* for Se 


X 2 for K e 


/i° L cal = 500 

/^° L cal = iooo 

/^° L cal = 5 000 


2.10 
3.02 x 10 
5.54 x 10 3 


2.12 x 10 
5.50 x 10 2 
1.31 x 10 6 



TABLE I: Results of the reduced x 2 test of the goodness of fit 
for Se and Ke calculated from the maps with different levels 
of non-Gaussianity as compared with the corresponding mean 
values obtained from MC simulated CMB input maps with 
0. 



A question that arises at this point is whether S and 
K indicators have sufficient sensitivity to detect devi- 
ation from non-Gaussianity for the values of the non- 
Gaussian parameter within the WMAP 7 yr bounds 
-fO < < 74 21]. Table |n] makes clear that for 

maps with /^L al within this interval, one has a negligi- 
ble value of x 2 , which makes it apparent that there is no 
significant overall departure of power spectra Si and Ke 
for /^L al = -10 and /jSjg* 1 = 74 from the correspond- 
ing MC mean power spectra obtained from the Gaussian 
( /nl = 0) ma P s J i n agreement with the nearly over- 
lapping symbols of Fig. @] This makes it apparent that 
the bispectrum based estimator that was employed by 
the WMAP team [t| [2l[ is more sensitive to primordial 
non-Gaussianity of local type than the S and K indi- 
cators used in the present paper (see also the related 
references [29]). Thus, for example, the deviation from 
Gaussianity as captured by our indicators S and K for 
simulated maps with /^L al = 74 is 4 orders of magnitude 
smaller than that for maps with f^ al = 1 000. 

To close this section, some words of clarification re- 
garding the forthcoming CMB data from the Planck mis- 
sion are in order. Planck combines high-angular resolu- 
tion and sensitivity with a wide frequency coverage (20 
GHz to 1 000 GHz) and will allow greatly improved fore- 
ground removal, thereby reducing many sources of non- 
Gaussian contaminants. The latest constraint on f^£ al 
from the WMAP-7 yr is = 32 ± 21 (68% confidence 

level), and from Planck we expect A/j!^ 1 ~ 5. This in- 
cludes cosmic variance and the detector noise [3(|. For 
a similar reduction at 2a confidence level, for example, 
the non-Gaussian procedure of our paper does not have 
enough sensibility to detect deviation from Gaussianity 
within the resulting Planck range of /^l • A prelimi- 
nary analysis shows that a new procedure with a different 
way of scanning input maps (through large pixels cells) 
seems to be capable of detecting such a tiny deviation 
from Gaussianity, i.e. within a narrower range of /^L al 
expected from the Planck mission. 



Non- Gaussian Parameter 


x" for Se 


X* for Ke 


Hocal i fi 

Jnl — — ru 
/N° L cal = +74 


2.80 x 10 -5 
1.10 x 10" 3 


4.80 x 10 -5 
1.00 x 10" 2 



TABLE II: Results of the x 2 test of the goodness of fit for 
Se and Ke calculated from the maps with /JjL al = — 10 and 
74 relative to the corresponding mean values obtained 



/local 
NL 

from MC simulated CMB input maps with /Jj 



= 0. 



/local 
NL 



B. Analyses and results for simulated and data 
maps 



Thus, TablclH along with Fig.[3l shows that the statis- 
tical indicators S and K can clearly detect deviation from 
Gaussianity of simulated maps endowed with primordial 
non-Gaussianity of local type with /^l 3,1 ^ 500. 



Besides the five frequency bands 7 yr maps released - 
K (22.8 GHz), Ka (33.0 GHz), Q (40.7 GHz), V (60.8 
GHz), and W (93.5 GHz)- the WMAP has also pro- 
duced a full-sky foreground-reduced ILC-7 yr map, which 
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FIG. 3: Low I average power spectra of skewness St (left panel) and kurtosis (right panel) Kt calculated from the Gaussian 
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500, 1 000, and 5 000) input simulated CMB maps. The la error bars are also indicated 



with a small horizontal shift to avoid overlap. The 95% confidence level relative to the Gaussian maps is indicated by the 
dotted line. 
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(/NL al = 0) and non-Gaussian simulated CMB maps for f^l* 1 = -10 and = 74. The 95% confidence level relative to the 

Gaussian maps is indicated by the dashed line. 



is formed from a weighted linear combination of these five 
frequency band maps in which the weights are chosen in 
order to minimize the galactic foreground contribution. 

The first-year ILC map has been explicitly stated as 
inappropriate for CMB scientific studies [28|], but in the 
subsequent ILC versions, including the ILC- 7 yr, a bias 
correction has been incorporated as part of the fore- 
ground cleaning process, and the WMAP team suggested 
that these maps are suitable for use in large angular scales 
(low £) analyses, although they have not performed non- 
Gaussian tests on these versions of the ILC maps pjl l3ll ? 
]• 

In a recent paper [l4| we have performed an analysis 
of Gaussianity of all the available full-sky foreground- 
reduced maps by using S and K non-Gaussianity indi- 
cators in the search for departure from Gaussianity on 
large angular scales 13]. We have shown that the full-sky 



foreground- reduced WMAP maps, including the ILC-7yr 
maps [l6[ , present a significant deviation from Gaussian- 
ity, in agreement with the results of Ref. [l3|, 0, H3 • 

An interesting question that arises here is how one 
compares the level of deviation from Gaussianity of the 
full-sky ILC maps with the level of primordial non- 
Gaussianity of local type with different amplitude pa- 
rameters /NL a1 - 8 We undertake this question by using 



We note that as the simulated CMB maps are full-sky maps, 
and sky cuts induce non-Gaussianity analyses. Thus, in order to 
make a comparative analysis of Gaussianity between simulated 
and data maps, a simple suitable choice of CMB data map is 
to take the equally full-sky data maps such as the ILC-7 yr. 
Nevertheless, in Appendix [A] we present a comparative analysis 
between maps with f^l al for the WMAP-7 yr bounds and the 
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simulated maps with the same smoothed 1° resolution of the ILC-7 yr WMAP map. The 95% confidence level relative to the 
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our indicators to carry out a comparative analysis of 
non-Gaussianity of simulated maps for the values of non- 
Gaussian parameters equal to the bounds reported by 
the WMAP team, and the ILC-7 yr map. 



Map 




X A for St 


X A for Ki 


rlocal 

Jnl 


= -10 


6.80 x 10" 2 


3.80 x 10" 1 


rlocal 

Jnl 


= +74 


7.10 x 10~ 2 


4.24 x 10 _1 


ILC 




8.96 


5.88 x 10 



TABLE III: Results of the reduced \ 2 test of the goodness 
of fit for St and Ki calculated from ILC-7 yr and MC sim- 
ulated maps for the 7 yr WMAP lower and upper bounds 
values of /jjl • in these calculations we have used simulated 
maps with the same smoothed 1° resolution of the ILC-7 yr 
WMAP map. These values quantify collectively the deviation 
from Gaussianity relative to MC Gaussian simulated maps 

(/NL Cal = 0). 

Figure [5] shows the low I mean power spectra of the 
skewness Si (left panel) and kurtosis Kg (right panel), 
calculated from simulated Gaussian (f^£ = 0) and non- 
Gaussian maps for /^L al = — 10, 74, which are the lower 
and upper bounds of the nonlinear parameter reported 
recently by WMAP [2l|. We have calculated the power 
spectra by using 1 000 maps for each value of /nl^- This 
figure also displays St and Kg computed from the full- 
sky foreground-reduced ILC-7 yr map. Figure O gives a 
clear indication that the deviation from Gaussianity of 
the ILC-7 yr map (detected by S and K) is greater than 
the non-Gaussianity of MC simulated maps for the lower 
and upper bounds of /-Jj£ al provided by WMAP [2l| . 



Q, V and W frequency maps. 



Clearly, in the calculation of the mean power spectra Si 
Ki of Fig. [3J we have used simulated maps with the same 
smoothed 1° resolution of the ILC-7 yr WMAP map. 

In order to compare quantitatively the deviation from 
Gaussianity, we have calculated the x 2 t es t °f t ne good- 
ness of fit for Si and Ki calculated from ILC-7 yr along 
with the simulated maps for /^l^ = —10,74 relative 
to the mean power spectra of Gaussian maps. We have 
gathered together in Table IIIII the obtained values for 
the reduced \ 2 - Figure [SJ along with Table ITTT1 shows, 
on the one hand, a significant deviation from Gaussian- 
ity in the full-sky foreground-reduced ILC-7 yr map, in 
agreement with the statistical analysis made in Ref. [lj] . 
On the other hand, it is apparent that the level of non- 
Gaussianity of the ILC-7 yr maps is higher than that 
of the simulated non-Gaussian maps with f^£ al equal to 
the bounds of the WMAP-7 yr data. This comparison 
provides indications of the suitability of the foreground- 
reduced ILC-7 yr map as Gaussian reconstructions of the 
CMB sky. 

V. CONCLUDING REMARKS 

The physics of the early universe can be probed by 
measurements of statistical properties of primordial fluc- 
tuations, which are the seeds for the temperature CMB 
anisotropies. Thus, the study of non-Gaussianity of 
these anisotropies offers a powerful approach to probe 
the physics of the primordial universe. It is essential, for 
example, to discriminate or even exclude classes of infla- 
tionary models, and also to test alternative scenarios of 
the primordial universe. 

Since one does not expect a single statistical estimator 
to be sensitive to all possible forms of non-Gaussianity, 
it seems important to test CMB data for deviations from 
Gaussianity by employing different statistical tools to 
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quantify or constrain the amount of any non-Gaussian 
signals in the data, and extract information concerning 
their potential origins. 

In recent papers [H, EH we proposed two new large- 
angle non-Gaussianity indicators, based on skewness and 
kurtosis of large-angle patches of the CMB sky sphere 
and used them to find significant large-angle devia- 
tion from Gaussianity in masked frequency bands and 
foreground-reduced CMB maps. 

Simulated CMB maps with an assigned primordial 
non-Gaussianity of a given type and amplitude are im- 
portant tools to study the sensitivity, power, and limita- 
tions of non-Gaussian estimators. They can also be used 
to calibrate non-Gaussian statistical indicators, and to 
study the effects of foregrounds and other non-Gaussian 
contaminants. 

In this paper we have addressed the question as 
to whether the non-Gaussian indicators proposed in 
Refs. [H, [HI have sufficient sensitivity to detect non- 
Gaussianity of local type, particularly with amplitude 
/^ £ al within the 7 yr WMAP bounds [H|. To this end, we 
have used our statistical indicators along with 6 000 simu- 
lated maps equipped with non-Gaussianity of local type 
with various amplitudes. From these simulated maps, 
which include the Gaussian one = 0), we have gen- 

erated 6 000 S maps and 6 000 K maps (see, e.g., Fig.[TJ), 
calculated the associated low £ mean power spectra Si 
and Ki, made a study of the sensitivity and strength, 
and determined the limitations of non-Gaussian estima- 
tors S and K. By using the mean power spectra of the 
simulated non-Gaussian maps along with the x 2 test of 
goodness, we have shown that S and K indicators can be 
used to detect deviation from Gaussianity of local type 
for typically /^ al > 500 (see Fig. [3] and Table HJ. Thus, 
our indicators do not have enough sensitivity to detect 
deviation from Gaussianity of local type with the non- 
linear parameter within the 7 yr WMAP bounds. This 
makes it apparent that the bispectrum based estimator 
employed by the WMAP team is more sensitive to pri- 
mordial non-Gaussianity of local type than the S and 
K. 

However, by using the same procedure, which is based 
upon the skeewness S and kurtosis indicators K of 
patches of the CMB sky sphere, we have shown that 
these indicators do not have enough sensitivity to capture 
non-Gaussianity when the nonlinear parameter f^ 31 lies 
within the 95% confidence level interval reported by the 
WMAP team [21] (cf. Fig. S and Table [ft]). The posi- 
tive outcome of the analysis performed in our previous 
works [l3|, [HI , together with the present outcome, seems 
to indicate that the deviation from Gaussianity captured 
in Refs.[l3l [l4| is not of primordial nature, although it 
might have a primordial component. 

Finally, we have also made a comparative study of non- 
Gaussianity of simulated maps and of the WMAP full-sky 
foreground- reduced 7 yr ILC map[l6|. which is summa- 
rized in Fig. [S] and Table Mil An interesting outcome of 
this analysis is that the level of non-Gaussianity of ILC- 



7 yr ( /nl ~ 775) i s higher than that of the simulated 
maps for /, 



local 
NL 



within observational bounds 211. This 



renders quantitative information about the suitability of 
the foreground-reduced ILC-7 yr map as a Gaussian re- 
construction of the CMB sky. 
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Appendix A 

Here we present the results of a comparative analysis 
of deviation from Gaussianity performed by using Q, V, 
W band maps and simulated maps equipped with non- 
Gaussianity of local type for f^ al = — 10 and 74, which 
are the WMAP- 7 yr bounds for fjffi? 1 . As the calculations 
are similar to those of Sec. IIV1 we refer the readers to that 



section for more details on how they are made. 

A word of clarification is in order before proceeding 
to the comparative analysis. In order to have an esti- 
mation of the role of the masking process in simulated 
maps with large values /Jf™ 1 , we have calculated the re- 
duced x 2 °f the mean Si and Ki for maps endowed with 
/k-cai = 500,1000,5 000 without and with the KQ75-7 
yr mask. These values quantify collectively the devia- 
tion from Gaussianity relative to MC Gaussian simulated 
maps (/^°£ al = 0). Tables US and H contain the results 
of our calculations. These tables show that the greater 
/^L al is the smaller the ratio is between x 2 values cal- 
culated for masked and unmasked simulated maps, re- 
spectively. Hence, the relative role of non-Gaussianity 
induced by the KQ75-7 yr mask in simulated maps is 
smaller for higher values of /J^" 1 , as we expected from 
the outset. Thus, the role of the masking process should 
be bigger for /^l 3,1 , taking the lower and upper bounds 
values of the WMAP-7 yr, an issue that will be discussed 
in the following. 

We begin the above-mentioned comparative analysis 
by recalling that one either has the full-sky contami- 
nated (Q, V, W) band maps to compare with full-sky 
simulated maps, or one masks both band and simulated 
maps. While in the former we have foreground maps 
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Full-sky map 


X* for S t 


X* for K e 


/iT 1 = 500 

tiT 1 = 1000 

fkT 1 = 5 000 


2.10 
3.02 x 10 
5.54 x 10 3 


2.12 x 10 
5.50 x 10 2 
1.31 x 10 6 



TABLE IV: Results of the reduced x 2 test of the goodness of 
fit for the mean power spectra Si and Ke as compared with 
the corresponding mean values obtained from MC simulated 
CMB input maps with /Jjf; al = 0. Full-sky simulated maps 
were used. 



Masked 


map 


X 2 for S e 


X 2 for K e 


/•local 

7nl — 


= 500 


3.40 x 10 


5.67 x 10 2 


/•local 

/nl — 


1000 


4.69 x 10 


6.81 x 10 2 


/•local 

/nl — 


5 000 


3.79 x 10 3 


4.79 x 10 5 



TABLE V: Results of the reduced \ 2 test of the goodness of 
fit for the mean power spectra Se and Ke as compared with 
the corresponding mean values obtained from MC simulated 
CMB input maps with /jj 3 ' = 0. The simulated input maps 
were masked with the KQ75-7 yr mask. 



Full-sky map 


X" for Se 


X 2 for Ke 


/•local i ri 

/nl — — iu 


2.80 x 10 -5 


4.80 x 10 -5 


/nT' = +74 


1.10 x 10~ 3 


1.00 x 10" 2 


Q 


8.79 x 10 14 


1.75 x 10 23 


V 


6.15 x 10 14 


3.17 x 10 23 


w 


3.14 x 10 13 


9.14 x 10 21 



TABLE VI: Results of the reduced x 2 test of the goodness of 
fit for St and Kt calculated from Q, V, and W band maps and 
MC simulated maps for the lower and upper bounds values of 
/JSl* 1 obtained from WMAP-7 yr data. These values quantify 
collectively the deviation from Gaussianity relative to MC 
Gaussian simulated maps (/nT 1 = 0). Full-sky simulated 
and frequency maps were used. 

that are quite contaminated, in the latter an induced 
non-Gaussianity arises from the masking process in sim- 
ulated maps. 

We have gathered together in Table lVll the values of the 
reduced x 2 for the full-sky Q , V, and W frequency maps 
along with the values of the reduced x 2 for maps endowed 



with non-Gaussianity for /^l = — 10 and 74. This table 
shows that the full-sky Q , V, and W band maps present a 
non-Gaussianity, as captured by our indicators, of several 
orders of magnitude higher than those of simulated maps 
for the lower and upper bounds values of f^ 3,1 obtained 
from WMAP-7 yr data. Clearly, this is not a surprising 
result since the full-sky Q, V, and W band maps are very 
foreground contaminated, but it shows the suitability of 
our skewness and kurtosis indicators to detect such a 
huge difference in the levels of non-Gaussianity. 

We have collected together in Table IVTTl the results of 
the calculations of the reduced x 2 for masked maps of the 



Masked 


map 


X 2 for St 


X 2 for K t 


.r local 
,/NL — 


-10 


0.79 


0.59 


/■local 

/nl — 


+74 


0.90 


0.85 


Q 




4.04 


3.20 


V 




0.86 


0.56 


w 




0.73 


0.60 



TABLE VII: Results of the reduced x 2 test of the goodness 
of fit for Se and Ke calculated from Q, V, and W band maps 
and MC simulated maps for WMAP-7 yr lower and upper 
bounds values of f^ 1 . The KQ75-7 yr mask was used for 
both frequency and simulated maps. These values quantify 
collectively the deviation from Gaussianity relative to MC 
Gaussian simulated maps (/nl = 0)- 



Q, V, and W bands, and for simulated input maps with 
/i^ al = -10 and 74. Tables Ell and [VTH show that, as 
expected, the KQ75-7 yr mask reduced significantly the 
levels of non-Gaussianity of bands maps, bringing them 
down several orders of magnitude. Even with the mask 
the detected level of non-Gaussianity for the Q band map 
is about 4 times the small level of the simulated masked 
maps for the WMAP-7 yr upper bound. As for the V 
and W bands and simulated maps, Table IVTTl shows that 
these bands maps present deviation from Gaussianity, as 
measured by our indicators, of similar order to those of 
maps with fj$Z ai in tlie WMAP-7 yr interval. 
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